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TECHNICAL MEMORANDUM X- 64920 


GRAVITATIONAL CLOCK: A PROPOSED EXPERIMENT 
FOR THE MEASUREMENT OF THE GRAVITATIONAL 

CONSTANT G 


I. INTRODUCTION 


Historically, the gravitational constant, G, was the first fundamental 
constant measured. Yet, today, it is the least accurately known of all funda- 
mental constants. This strange fact, until now, has been the result of two 
reinforcing circumstances: (l) the gravitational force is relatively very small, 
and (2) the gravitational constant is not necessary (important) in most prac- 
tical physical problems; e. g. , in electronics, or in planetary physics where G 
is not needed, and GM which can be measured astronomically very accurately is 
used. The irony of this situation is that since gravitational forces are small 
and purely mechanical, measurement of G is possible, in principle, to even 
greater accuracy than those constants that graze the domain of the quantum. 
Nevertheless, experimental obstacles have kept the accuracy of G to 1 part in 
10 3 . However, through the use of space-age technology and satellite-borne 
experiments, most of these difficulties may be overcome, but this is not 
enough. All this would not be important except that the physical importance of 
knowing G very accurately has grown significantly in the last few years. The 
impetus for this knowledge is twofold: (l) the desire to test competing gravita- 
tional theories, and (2) experimental evidence of Van Flandern that G is 
changing with time [ l] . 

Based on new and old observational data, Van Flandern concludes that 
from anomalous secular acceleration of the moon’s mean longitude, the 
fractional rate change of G is 


^ = _(8 ± 5) x 10 _11 /year 
G 


(1) 



Two alternate, but much less probable, explanations are a uniform expansion 
of all space or a solar system mass loss (presumably from the sun). The 
latter possibility is not consistent with actual solar mass losses, whereas a 
uniform expansion is ruled out (tentatively) since this implies a continuous 
increase in angular momentum; e.g. , of the earth which is also not observed. 
Thus, the manifestly non-Einsteinian general relativistic result that G =£ 0 
remains. However, it is a character of the Brans-Dicke-Jordan scalar-tensor 
theories that G can be nonzero [2,3] . The same is true in the Vector-metric 
[4] or two-tensor metric theories [5] in which G is related to the second field. 
However^ - Van Flandern indicates [l] that his value of G [equation ( l) ] - is more 
closely related to the cosmologies of Hoyle and Narlika [6] or that of Dirac 
[7] than with the Brans-Dicke theory. Van Flandern’ s contention about Brans- 
Dicke theories does not seem to be correct, since the value of G predicted by 
any scalar-tensor theory is strongly dependent on the value of poorly measured 
cosmological parameters such as the universe mass density and the decelera- 
tion parameter [8] . We must also contend with the possibility that there may 
be spatial variations in G such as occur, for example, in the scalar-tensor 
theories [9] or in Ni’s new theory [10] . In general, this second type of varia- 
tion of G can be summarized using the parameterized post-Newtonian (PPN) 
formalism [11] as indications of preferred frame effects or possible noncon- 
servation of total four-momentum [9] . Preferred frame effects seem to be 
ruled out on the basis of experiment [9]. However, the experimental limits on 
some aspects of nonconservative theories in the PPN approximation have not 
been established, since these effects seem to be at or beyond our present 
experimental capabilities [9]. 

In the new PPN parameters, the spatial variation of G can be written as 


G = 1 - (4/3 -y-3 -£ 2 )U 


( 2 ) 


where U is the gravitational potential, j 3 and y are the classical Robertson 
parameters (both equal to one in Einstein general relativity) , and t, 2 is one °f 
four PPN parameters that indicates a possible nonconservative theory if non- 
zero [ 9] . In general, t, 2 indicates a breakdown of the equivalence principle 

(M . =£ M ), [Possible preferred frame effects have been neglected 

passive active & 

in equation (2) . ] In the Brans-Dicke theory [9, 12] , /3 = 1, y = (03 + l) /(o> + 2) , 

where co is the usual Brans-Dicke coupling parameter and — 0 so that 


(4)3 - y _ 3 - ^ 2 ) = l/(w + 2) 


(3) 


2 



which then leads to spatial variations in G. However, these variations are still 
three or four orders of magnitude smaller than the best experimental accuracy 
of G [13, 14]. Recently Smalley 1 has shown that nonconservative theories, in 

particular theories for which ^ , v * 0 > are possible; i.e. , consistent with 

the field equations (Bianchi identities), gauge conditions on the metric, and the 
Newtonian limit of the conservation laws. One such theory, which is a modifica- 
tion of the Brans -Dicke theory, predicts the same PPN parameters with the 
exception that 


r = so 

(2u) + 3) a ’ 

where a is a coupling constant in the expression for the nonzero divergence 
of the energy-momentum tensor 


<b 

V'JA‘ 87T -JJ. \V 


(5) 


and <p is the scalar field. This theory agrees in all aspects, except for t, 2 + 0, 
with the usual Brans-Dicke theory. It should be contrasted with the original 
nonconservative scalar-tensor theory of Jordan [15] but does not, however, 
suffer from the criticisms of Fierz [16] or Bondi [ 17] , at least at the level of 
PPN approximation. By comparing equations (3) and (4), it is possible that in 
a nonconservative Brans-Dicke theory the spatial variations in G could dis- 
appear; i.e. , for the right choice of cr, G = 1 in equation (2). Thus, the viola- 
tion of the strong principle of equivalence can be masked if the theory is also 
nonconservative, at least on the level of the PPN approximation. Thus we note 
that there are two competing effects which lead to variations in G — spatial and 
temporal variations — and, according to our discussion above, it is not certain 
whether a clear separation of the effects can be made. 

Suppose we have a detector in orbit about the earth and we attempt to 
measure G. Since the gravitational force is small, we know, a priori, that the 
experiment will be long term, perhaps lasting for months or a year* However, 


1* Smalley, L* L*: Gravitational Theories with Nonzero Divergence of the 
Energy -Momentum-Tensor* To be published in The Physical Review. 
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the eccentricity of the earth's orbit about the sun will produce a 1 percent 
modulation of the solar gravitational potential U at the position of the earth, 
which could then produce a yearly spatial amplitude variation in G of 2 


AG 


G 


3 x IQ- 11 


( 6 ) 


This variation is of the same order of magnitude as the temporal variation of 
G reported by Van Flandern. Therefore, any experiment that can determine 
G at the level of accuracy of 1 part in 10 11 must contend carefully with these 
competing effects. One possibility is a solar satellite with nearly zero eccen- 
tricity so that the temporal variations dominate or vice versa. 

Finally, one last effect must be discussed. It has been tacitly assumed 
in this discussion that the Newtonian force law is an inverse square law. 
Although this may appear to be the same problem, it is not, but depends on the 
level at which general relativistic corrections become important. However, 
this test suffers identical experimental problems of accuracies of 1 part in 10 4 
on the earth primarily because of errors in location of test masses [ 14, 18] . 
Until recently it was thought that the best test of the inverse square law came 
from the astronomical observations of the advance of the perihelion of Mercury. 
Long and Ogden have shown [ 19] that by using a double-ring source, very uni- 
form gravitational fields can be produced in the laboratory. This means that if 

they consider the historic inverse-square-law failure in the form r - ^ + ^ , the 
advance of the perihelion of Mercury yields the limit, <j> < 4 x io -9 , whereas the 
use of their uniform laboratory fields could conceivably give a limit of M 3x 
10- 12 . This stronger limit is, of course, directly related to the relatively 
lowered importance of positional accuracies for the test masses because of the 
highly uniform laboratory field. 

Thus, it seems that an inverse-square-law failure could mask a spatial 
variation in G. Experimentally three intertwined effects remain that could con- 
ceivably introduce effects at the same order of magnitude. However, these 
three effects combine to enhance the theoretical importance of measuring G to 
higher accuracy if we hope to establish the correct gravitational theory. Most 
likely, these are questions we must answer before we can solve the important 
problems of initial data, quantum gravity, or unified field theory. 


2. Weiss, R.: MIT internal report. 
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In Section II, recent attempts or proposed experiments for the deter- 
mination of G are reviewed. However, recent work on null-experiments is 
omitted; i.e. , Eotvos experiments that test the equivalence principle, except 
experiments that could conceivably be used to measure G directly. A satellite- 
borne gravitational clock experiment from which G can be obtained directly is 
proposed in Section III, and conclusions and some of the prospects for direct 
measurement of G are described in Section IV. 


II. PRESENT STATE OF THE ART FOR OBSERVATIONS OF G 


It is one thing to say that G has been measured to 1 part in 10 10 , and it 
is entirely another to say that G/G is known to the same value. The latter does 
not imply the former, not only because the value is model dependent, as in Van 
Flandern's work, but because a known variation in G never implies a knowledge 
of the magnitude of G. This becomes even more apparent when we recall that 
G is known at best to 1 part in IQ 2 * 4 after several centuries of experimental and 
theoretical research [13]. 

Basically there are two ways to measure gravitational forces (following 
the outline of Weiss, previously cited in footnote 2) : (l) by use of static 
mechanical or electrical forces such as springs, fibers, or electromagnetic 
devices; or (2) by gravitational inertia reaction forces. In the first category 
are the classical torsion balances of Cavendish and Eotvos and their modern 
improvements [20,21] as well as assorted electromechanical devices such as 
gravimeters, accelerometers, and gradiometers. The latter type of experi- 
ments can be generally classified as free-fall experiments, with the exception 
of the Beams' balance experiment [14,22] or the proposed rotating flat-plate 
experiment of Berman [23] . The primary difficulty with inertia experiments 
on the earth is suspension in the earth's gravitational field, which tends to 
overwhelm the dynamical effects. 

Examples of free-fall inertia reaction experiments that have been pro- 
posed to measure G are: 

1. Measuring the period of synchronously orbiting spheres by Berman 
and Forward [24] or Chapman (from the work previously cited in footnote 2) . 

2. Measuring the period of oscillation of a test mass through a bored 

hole in a larger spherical mass by Berman and Forward (from the work 

previously cited in footnote 2) and Berman [ 25] , by a modification of the 
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proposed Worden and Everitt 3 Eotvos experiment, or by the flat-plate spherical 
mass oscillator proposed here. 

3, Tangential acceleration balance by Blood (from the work previously 
cited in footnote 2), which is essentially a modification of the Beams’ balance 

method. Weiss (see footnote 2) proposes to use the same experimental 

arrangement in a ’’clock” mode. 

4. Centrifugal balance systems by Wilk (from the work previously 
cited in footnote 2) , or through a modification of the rotating Eotvos force- 
balance proposed by Chapman [26,27] . 

We will touch, without mathematical details, on most of these experiments, 
stating briefly their methods, advantages, and disadvantages as well as their 
ultimate expectations. 

As far as this author knows, there appear to be at present only two 
ongoing experiments to measure G: the ordinary, but supposedly improved 
sensitivity and experimentally isolated, torsion experiment of Cook and 
Marussi [20,21] and the constant angular acceleration experiment of Beams 
[ 14, 22] . Also included here should be the active proposal of Worden and 
Everitt ( see footnote 3) on which work on subsystem instrumentation is pro- 
ceeding. However, this experiment is designed as an E6tvos experiment and 
would require modification along the lines of Weiss’ modification, discussed 
below, of Blood’s modified Beams' balance experiment (see footnote 2). 

The standard torsion experiment (e.g. , see Reference 18) seems to be 
limited by mechanical noise to accuracies somewhat less than 1 part in 10 6 , 
which Cook quotes as the limitation of their measurements of masses, lengths, 
and time. The Beams inertia reaction experiment (Fig. l) features a rotating 
table accelerating with a constant a ng ular acceleration. On the table are two 
large masses in a dumbbell arrangement which, when rotated about an axis 
perpendicular to their line of centers, produce a torque ( T) on a small cylin- 
drical bar. The experimental procedure is to compare the measured torque 
produced on the small cylinder with and without the large masses on the 
angularly accelerated table. The difference between these two torques allows 
one to calculate G directly. In their initial experiment, the small cylindrical 
mass was suspended by a quartz fiber, and ordinary bearings were used for 
the accelerating table [ 14, 22] . Thus far, mechanical noise has limited the 
experiment to accuracies of about 1 part in 10 4 * . The limiting accuracy seems 


3. Worden, P. W. ; and Everitt, C. W. F.: Test of the Equivalence of 

Gravitational and Inertial Mass Based on Cryogenic Techniques. Preprint, 

Stanford University, Jan. 1973. 
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Figure 1. Beams’ balance configuration (two large masses are constrained 
to accelerate with angular acceleration, a. , about a small cylinder; from 
the torque, T, sensed through the twisting of a quartz fiber about the 
vertical, the gravitational constant, G, can be calculated). 


to be related to the certainty in the 
location of the center of mass of the 
two large spheres of about 5 parts 
in 10 s , This then results in an ulti- 
mate design capability of 1 part in 
10 G , which they believe they can 
achieve providing they use magnetic 
suspension systems and/or air 
bearings. 

A novel example of a pro- 
posed experiment for determining G 
is the rotating flat plate [23l de- 
picted in Figure 2. In this experi- 
ment two optically flat and parallel 
rectangular solids are arranged so 
that one is suspended from a quartz 
fiber and is, therefore, free to 
oscillate with amplitude 8 about the 
vertical. The other plate is then 
driven with a uniform angular 
velocity about the same axis. The 
suspended plate then resonates 



Figure 2. An optically flat plate rotates 
below another flat plate that is sus- 
pended by a quartz fiber (the gravi- 
tational constant can be obtained 
from the period of oscillation 
of the suspended plate about 
the vertical) . 
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with some period from which G can be mathematically extracted providing 
anharmonics can be eliminated. By careful design, the fourth-order terms in 

0 can be removed from the potential, but the forces (i.e. , torques) are very 
small; and since the excitation mechanism depends on the second harmonic of 
the gravitational field, the sixth-order terms are important. Although physical 
dimensions can be fixed extremely accurately because of the nice optical prop- 
erties of the plates, and since the system would be highly evacuated, it is 

— purely a mechanical system and, therefore, mechanical noise limited. Thus, 
because the forces are so small, the accuracy of the system is limited by 
anharmonics and noise. Berman [23] estimates ultimate accuracies of about 

1 part in 10 5 . Although not discussed by Berman, it is obvious that this 
experiment could operate in the constant angular acceleration mode similar to 
that of the Beams’ balance. We have not as yet analyzed this possibility, but 
its simplicity is appealing. 

In the first category of inertia reaction experiments are the miniature 
orbital systems such as, for example, that proposed by Berman and Forward 
[24] of two identical spherical masses, M, of radius r, in close orbit a dis- 
tance R i. 2r apart as shown in Figure 3. A realistic period that depends on 
the density of the spheres and R is about 100 minutes. Although this system 
allows one to calculate G directly from the elementary centripetal equation 


S2 2 = 


2GM 

~w 


(7) 


by measuring the orbital frequency, , the orbital radius, R, would be poorly 
determined and would be a major source of error limiting precision to about 
that of the older mechanical systems. Serious sources of error and limitation, 
not discussed for this particular experiment by Berman and Forward, are 
gravity and centrifugal gradients (with their own periods) that will tend to 
distort, nontrivially, the orbital motion of the spheres and their periods of 
rotation depending on the orientation of the orbital plane of the spheres with 
respect to the orbital plane of the satellite housing the experiment. We 
conclude, as they do, that this is not yet a practical laboratory experiments 

The second category of inertia experiments is the standard physics 
textbook problem, shown in Figure 4, of a test mass in free fall through a 
tunnel bored in a much larger sphere of uniform density. Such a free-fall 
system has also been investigated in detail by Berman and Forward [ 24] , and 
there is even a modification by Berman [25] , shown in Figure 5, that could be 
done on the surface of the earth by using a levered beam balance. However, 
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Figure 3. Two spheres of mass M and radius r mutually 
orbit each other at a radius R 2r. 



FALL 

MASS 


Figure 4. A small test mass oscillates in a tunnel bored through 
a larger uniform density sphere. 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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Berman concludes that the earth- 
based experiment would not improve 
on the accuracy of G over the usual 
mechanical systems primarily 
because of gravity gradients and very 
serious alignment problems that 
would alter the period of oscillation 
significantly. However, the free-fall 
system of a test mass oscillating in a 
tunnel bored through a larger sphere 
bears more promise. 

Again, this gives a direct 
measure of G from the period of 
oscillation of the test mass. How- 
ever, because of nonspherical 
symmetry ( after the tunnel is bored) , 
the mathematics is no longer the 
simple theoretical problem originally 
proposed. In fact, the once obvious 
advantage of spherical symmetry is 
lost. Another limitation is the 
requirement for large amplitudes ; 
otherwise there is no physical need 
in the first place for a tunnel 
through a large sphere. However, 
large amplitudes exacerbate gravity 
gradients. The experimenters thus conclude that the only allowable orientation 
for the experiment is perpendicular to the orbital plane of the satellite along 
which axis the gradient forces are constant. Even so, the accuracy of the 
system is ultimately limited by determination of the radius and density 
variation of the larger sphere to about 1 part in 10 6 . 

Worden and Everitt (see footnote 3) have proposed a free-fall Eotvos 
experiment in which a small cylinder is centered and parallel inside a larger 
cylinder, as shown in Figure 6. The acceleration which tends to decenter the 
masses would be a measure, among other possibilities, of a violation of the 
equivalence principle. Experimentally, no violation has been observed [28,29] 
down to an Eotvos ratio of 1 part in 10 11 . However, their system, which is to 
be cryogenic and magnetically suspended and centered by superconducting 
coils, is theoretically capable of detecting an Eotvos ratio of 1 part in 10 19 but 
is experimentally limited at present by helium management to about 1 part in 



Figure 5. An earth-based version of 
the test mass oscillating in a tunnel 
of a larger sphere (the acceleration 
of gravity, g, is compensated for 
by a fulcrumed balance so that two 
small test masses can oscillate 
inside tunnels bored through 
the larger sphere) . 
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Figure 6. A small cylindrical test mass is free to oscillate along 
the axis of a larger but hollo-w cylinder. 

10 14 , although they hope to improve this significantly in the near future by a 
better design of this subsystem. It is important to note that their system 
could be utilized in the oscillatory mode instead of using the acceleration sig- 
nal in a feedback loop to a magnetic controller that continually centers the tvvo 
masses. In fact, the controller system is a source of error since it introduces 
an oscillation about the null which then constitutes the minimum detectable 
reading. Fortunately, this error will be two or three orders of magnitude 
less than thermal noise. All other sources of error — such as gravity gradi- 
ents, orbital eccentricity, differences in homogeneity of masses that cause a 
noncentering of the geometrical center with the center of mass (CM) , non- 
coincidence of CM with the center of gravity (CG) , asymmetrical satellite 
masses, residual drag-free satellite oscillations, magnetic forces, and 
residual gases — , they claim, can be detected and controlled so that thermal 
noise becomes the limiting source of error. However, an orbital eccentricity 
does not produce a driving harmonic oscillation of their system since such 
accelerations are proportional to amplitude in the same way that gradient 
forces are proportional to amplitude. Thus, looking at the Fourier component 
of the acceleration at orbital frequency as a source of error as in an Ebtvos 
experiment will not work for the clock mode of operations. It is shown in the 
next section, in a discussion of the flat-plate spherical oscillator, that the 
system is essentially equivalent to a periodic Mathieus’ equation instead of a 
driven harmonic oscillator. However, this is only a mathematical detail. 
Nevertheless, their experiment is well developed and has the immediate advan- 
tage of being able to perform two experiments, although the oscillator experi- 
ment has geometrical limitations which will be discussed later. 
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The third type of free-fall inertia experiment is a proposal by Weiss 
to modify Blood's tangential acceleration balance (see footnote 2) to a clock 
mode instead of uniformly accelerating one set of masses with respect to the 
other set (as in the Beams' balance). The tangential balance clock is basically 
a double cruciform structure (Fig. 7) with the central three masses, called 
the rotor, oscillating through tunnels bored in three larger exterior stator 
masses. Amplitude anharmonics can be reduced by design below the detection 
threshoTd for a proposed laser interferometer detection system capable of 
measuring fractional periods to 10 -8 . Although thermal noise is considered the 
fundamental limit to precision, Weiss claims that suspension problems limit 
the experiment to an ultimate measurement of G to 1 part in 10 8 . 

The fourth and final category of inertia reaction experiments is the 
rotating Ebtvos centrifugal force-balance system of Wilk (from the work 
previously cited in footnote 2) or Chapman [26,27], which could also be 
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Figure 7. In the double cruciform structure the rotor masses oscillate 
through tunnels bored in stator masses. 



operated in the clock mode for direct determination of G. In this type of 
satellite Eotvos experiment, a small proof mass is constrained to oscillate 
about the null position either mechanically or electromagnetically, as shown 
in Figure 8. In the clock mode this constraint is removed so that the small 
mass is free to oscillate along the accelerometer’s sensitive axis under the 
mutual gravitational force between it and the larger mass. The right angular 
velocity of rotation of the larger disc mass and gravitational forces between 
the masses and cyclic gravity gradients then combine to restore the proof 
mass periodically toward the null position. The immediate advantage of this 
rotating system is that it attempts to compensate for gravity gradients through 


* LOCAL 
j VERTICAL 

i 

i ACCELEROMETER 



Figure 8. The proof mass oscillates along an axis of a slot cut from a 
larger disc mass which itself rotates about an axis perpendicular 
to the orbital plane that is in the plane of the paper. 
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centrifugal forces. The system is then described by a Mathieus 1 equation [30l 
of the form 


d 2 x 

w 


+ 


(a - 2q cos 2£)x = 


0 


( 8 ) 


where a & 0 , and q £ 0 ( see Section m) . Unfortunately, its solution is only 
marginally stable since for small but negative a and q £ 0 , it is near the 
border of the first instability region for the Mathieus 1 equation. If the system 
can be successfully operated in a stability region, the system is inherently as 
accurate as that of the double cylinder system of Worden and Everitt, provided 
cryogenic techniques and electromagnetic suspension are used. It should be 
understood that one does not neglect gravity gradient in the double cylinder 
experiment but that the system is easily operated in a stable mode, as will 
become obvious in the discussion of the flat-plate spherical mass oscillator in 
Section in. Another disadvantage of the rotating system is the method by 
which the system is rotated. This requires that the system be attached to a 
very stable, directionally oriented platform, whereas the cylindrical system 
or the flat-plate system is automatically stabilized along its own inertial axis 
and, thus, can fly ’’free” inside a larger "drag-free" satellite. 


III. FLAT-PLATE SPHERICAL MASS OSCILLATOR 


Consider a cylindrically symmetric flat plate of radius c and height 
2t with a circular hole of radius b < c bored along an axis perpendicular to 
and centered on the flat ends of the cylinder. We then allow a spherically 
symmetric mass of radius b' <; b to oscillate along this axis as shown in 
Figure 9. With the exception of the orbital mass systems, the combination of 
cylindrical and spherical symmetry is the simplest mathematical configuration 
that includes the effects of gravity and centrifugal gradients which can be used 
experimentally. Unfortunately, the more highly symmetric orbital system 
has severe experimental problems in determining the orbital radius, so that 
we are forced to consider a less symmetric system. Even the tunneled sphere 
problem depicted in Figure 4 is less symmetric because the tunnel destroys 
the useful spherical symmetry of the larger spherical mass. The cylindrical 
oscillator of Figure 6 has yet a lower symmetry, which implies a more 
untractable mathematical problem. However, these last two systems suffer 
experimental limitations that are somewhat circumvented by the inherent 
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geometry of the orbital, flat-plate, 
spherical mass system. In the 
case of the totally cylindrical 
oscillator (Fig. 6) , the gravita- 
tional restoring force tends to be 
smaller. This can be understood 
by noting that for a small displace- 
ment 6 x from the equilibrium 
position shown in Figure 10, the 
restoring force comes primarily 
from a slab of thickness 25 x at 
the opposite end of the outer cylinder 
which then produces a relatively 
weaker restoring force. In the 
(ideal) large spherical system, the 
restoring force, as is well known, 
is due to the mass contained inside 
the mathematical sphere whose 
radius equals the test mass dis- 
placement from the center of the 
sphere along the tunnel. One, thus, 
has large excursions, which in- 
creases the relative importance of 
centrifugal and gravity gradients 
that increase with the radius and 
also are modulated by a time- 
dependent periodic function at 
twice the orbital frequency for circular orbits. If the symmetry were exactly 
spherical, this would cause no problems, provided a stable mode is possible, 
since then the gravitational restoring force would also increase with radius. 
However, a tunnel of finite radius produces anharmonic amplitude effects 
which then require small amplitudes, and this implies a smaller attractive 
mass than for the flat-plate spherical mass oscillator. That is, for the flat- 
plate for an amplitude Z < 1, the attractive mass decreases as z , whereas 
for the sphere, it decreases as z 3 . Thus, there appear to be both theoretical 
and experimental advantages for the flat-plate spherical mass oscillator over 
the other two oscillators described. However, this statement must be qualified 
to the point that the amplitude of oscillation remains small, but this results in 
another advantage. 



Figure 9. Side view of a flat-plate, 
spherical mass oscillator (the 
sphere of radius b r oscillates 
through a tunnel of radius b 
bored through a tunnel in the 
flat-plate cylinder of 
radius c) . 
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Figure 10. For a displacement of 5x of the small interior cylinder, the 
restoring force comes from a small slab of the outer cylinder of 
thickness 26x a distance L-5x away from the center of 
mass of the small cylinder. 

To demonstrate the mathematical simplicity of the flat-plate spherical 
mass oscillator, consider first Figure 11 which depicts a spherically sym- 
metric mass of uniform density and radius r constrained to move along the 
axis of symmetry (z axis) of a thin hoop of radius p and surface density X. 
The mass elements for the hoop are then 


dM, = Xpdpda 
h 


(9) 


where a. is the angle of the differential of mass, dM^ , around the hoop, and 
the differential of mass for the sphere is given by ‘ 


dM = ix r 2 dr sin 0d Qd<p 
s 


( 10 ) 
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Figure 11. Geometry for calculating the gravitational potential between a 
thin hoop and a spherical mass on the axis of the hoop. 

where 9 is the azimuthal angle measured from the line of length 5, and <f> 
is the angle of rotation about this axis. We have then 


5 2 = p 2 + z 2 , (u) 


where z is the distance between the centers of mass and the distance D 
between mass elements is described by 


D 2 = 6 2 + r 2 - 26 r cos 9 


( 12 ) 


The differential gravitational potential experienced by dM 

h 

d <fr = -dM G f ^ Sin — 9d ^ dr = GM s dM h 
k sphere D ® 


is then given by 

(13) 
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That is, we can treat the sphere as a point particle of mass as long as 

the sphere is external to the hoop or to the generalization of the hoop into a 
short, hollow, right circular cylinder (Fig. 12) of height 2 t , centered at 
z = 0 , and which has external radius c and internal radius b > b' where b’ 
is the radius of the sphere. In this case, the potential at z becomes 



2-nv\.pdpd£ 

6 


* 


(14) 


where X is now a uniform volume density and £ is the variable of integration 
along the z direction and 


6 2 = (z - 4) 2 + P 2 


(15) 


After the integration over p , we obtain 


<A(z) = 2TTAGM 

s 


z - T / 

f yJ c 2 + x 2 

Z + T 



(16) 



Figure 12. Geometry for calculating the gravitational potential between 
a cylinder and a point mass on the axis of the cylinder. 
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The force on 


M 

s 


is then given by 


F s<‘> - -** . - -t 


= 2ff\GM g c 2 + (z - t ) 2 - N/b 2 + (z - t ) 2 

_ *Vc 2 + (z + t ) 2 + ^ b 2 + (z + t) 2 J 


« 27t\GM tJ 

s ' 


\4 b 2 + T 2 ‘/c 2 


=7) 


b 2 


(c 777 )% (b 2 + r 2 ) y 2 


(17) 


where we have expanded for small z. To compare the effect of the z 3 term, we 
consider the case where 10 2 r = 10b = c which yields the ratio of the third-order 
term to the second-order term. 


~ 5 x 10' 3 z 2 


(18) 


Thus, to reduce anharmonics (in this example) to an effect on the force of less 
than 1 part in 10 11 , we must keep amplitudes z £ 10“ 4 cm. Although this 
amplitude is small, capacitive or magnetic -inductive tracking (see footnote 3) 
appears sensitive and accurate to 0. 05 A so that small amplitude oscillations, 
we will find, permit a convenient mode of operation. Thus, we propose to 
operate the clock with the amplitudes that allow us to neglect the third-order 
term. 


Define the quantity f , which is the ratio in equation ( 17 ) of the first- 
order term to z , 


f 

c 


4x\GM T 
s 



1 


J7T7 


) 


(19) 
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In terms of the mass of the plate 


M 

P 


27t\t (c 2 - b 2 ) 


then 


f 

c 


2GM M 
P s 


(c 2 rwj 




( 20 ) 


( 21 ) 


In the two-body central force problem, 

mass p = M M /( M + M ) so that 
p s p s 

f 2G(M + M ) , 

.2 _ C _ P s' / 

" h=‘^= 


For large c and X = p = 16 g/cm 3 , w 2 ** 4^pG ; this implies a period of 
approximately 29 minutes. 

Consider now the same oscillator in orbit with its axis of oscillation in 
the plane of the orbit. The total force on the two masses is due to three 
things: (l) the force due to each other, (2) the force due to the earth, and (3) 
gradient force due to centrifugal forces. This latter force, although a fictitious 
force, must be considered since we are in actuality dealing with an orbital 
problem for which the centripetal acceleration is given by 


the frequency is related to the reduced 


( oo\ 


a = Rft 2 
c 


(23) 


where R is orbital radius and ft is orbital angular velocity. This implies 
that for small separation z. from the center of mass, the ith particle of mass 

m. will experience a repulsive centrifugal gradient force. 


F. = m z ft 2 

l i i 


( 24 ) 
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in the plane of the orbit. In the case of circular orbits 


J2 2 = 


GM 

IF" 


( 25 ) 


where M is the mass of the earth. The total force on each part of the 
oscillator becomes 


GM M 

M 1 z 1 = f c (z 2 - z x ) + -^-cosa + 


GM M 

M 2 " 2 = f c (z i - Z 2 } + “S*“ COsS3 + M 2 Z 2 n2 


(26) 


where we have used the geometrical quantities defined in Figure 13, and the 

subscripts 1 and 2 refer, respectively, to the flat-plate and the spherical 

mass. In terms of the central force problem defined by the coordinate 

z = z - z and the reduced mass m = M M /(M + M ) , we find that the 
1 2 J- 2 -L 2 

central force problem can be described by the single equation 


p z 


-f z + GpM 
c 


cos a cos /3 \ o2 

-5T -- sr ) + 

1 2 / 


(27) 


Expanding R and R in terms of z /R and z /R, respectively, we obtain, 
12 12 

after dropping higher order terms. 


* 3GjxM 2 

p.z = -f z + — z cos' 1 ilt 
c R 


L 3GM^\ , 3GMM _ n . 

VliH " + IS 3- 


(28) 
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M (EARTH) 


ORBITAL PLANE 



Figure 13. Geometry for calculating the central force equation between two 
masses, M and M , in the presence of an earth of mass M with orbital 

1 a 

angular velocity ft (all three masses lie in the orbital plane). 
Letting ft t = jj , we then obtain the equation 



t \ 


(29) 


which can be compared with the general Mathieus 1 equation given in equation 
(8) with 


a 



q = 


3 

4 


( 30 ) 
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Because the driving force in a Mathieus 1 oscillator is itself proportional to 
amplitude, it is not stable for arbitrary values of a and q as is the ordinary 
driven harmonic oscillator. For q=3/4, stability requires that 


-0.27 < a < 0.18 


(31) 


for the first region of stability. The second region of stability begins at about 
a > 1.72 which represents an upper bound (as if the boundary could be repre- 
sented by a straight line) . From equation (30) , we find 


1.23 



1.68 


w* 

0 * 


> 3.22 


(32) 


and so forth. For a given set of masses (and geometry) , the stability regions 
depend on the orbital radius through S2 ; these regions then oscillate with 
altitude. We now note the unusual consequence of the Matirieus’ equation: 
stable operation requires that the oscillator frequency be greater than the 
orbital frequency; i.e. , the character of the solution is very much different 
from that of a driven harmonic oscillator. However, within the stable ranges, 
it is possible to extract G from knowledge of both the period and the amplitude 
of oscillation. 

Another mode of operation occurs when the oscillator axis is perpen- 
dicular to the orbital plane as depicted in Figure 14. In this case 



(33) 


is an ordinary harmonic oscillator with frequency given by 


cor 



CO 2 + fi 2 


(34) 
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Figure 14. The masses oscillate perpendicular to the orbital plane 
as opposed to being in the plane as in Figure 13. 

Since the quantity S2 can be accurately determined from orbital tracking data, 
G can be easily extracted from equation (34) after measuring co^ . On the 

other hand, if the experiment is operated at several different altitudes, 
knowledge of (R) allows one to extract both G and M from the data. 

Thus, if we include this mode of operation along with that of the oscillator axis 
in the orbital plane, there are three independent ways to extract G from the 
data. 


It should be pointed out that the above calculations assume perfectly 
circular orbits. Noncircular orbits introduce orbital frequency harmonic and 
subharmonic terms into the force described by equation (28). The largest 
perturbation is at orbital frequency instead of the twice orbital effect given in 
equation (28) but is proportional to the orbital eccentricity, e . Thus, for 
circular orbits with e < 10 -8 , we can neglect these effects on frequency and 
amplitude measurements for the oscillator. 

Because of geometrical considerations, the radius of the sphere will 
not be known as accurately as the dimension of the flat plate, which can be 
fabricated more accurately. However, the uncertainty in the radius of the 
sphere (i.e. , the mass of the sphere) because of spherical symmetry can be 
transferred into an uncertainty in z. For a small uncertainty of about 1 part 
in 10 6 , which seems feasible, this contribution to the period is negligible. In 
principle, the period is simply related to the total masses and the overall 
geometry. Nonuniform density does, however, produce an uncertainty in the 
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position of the CM of the sphere and flat plate which tends to produce a torque 
because of the nonuniformity of the gravitational potential of the flat plate over 
the volume of the sphere* This torque then produces a periodic rotation of the 
sphere with respect to the plate which is a direct determination of the inhomo- 
geneity of the sphere. Very accurate determination of these rotations is 
possible through the detection of the London moment produced by a rotating 
superconductor. 

To monitor the oscillations of the sphere, we propose to enclose the 
spherical mass inside a thin spherical shell rigidly attached to the flat plate. 
Because of the spherical symmetry of the shell, the shell does not affect the 
motion of the spherical mass. The oscillator will then be constrained to move 
along one dimension by using a capacitance controller attached to the shell 
that can also record the amplitude and frequency of oscillations along the active 
axis. Superconducting loops can also then be attached to the shell, and by use 
of Josephson bridges, the rotation of the sphere can be measured. 

The measuring system introduces so-called tracking errors which, for 
instance, in an Eotvos null system would produce a feedback oscillation that 
would then represent the minimum detectable motion of the system. Worden 
and Everitt (see footnote 3) have discussed this iype of error and conclude 
that it is much less than amplitudes expected from thermal noise which can be 
approximated by 


<Ax> 



(35) 


where k is the Boltzmann constant and T is temperature. For a period of 
1 hour, M = 10 kg, T = 3k, then <Ax> ~ 10 -7 cm. The minuteness of the 
thermal errors compared with a signal of 10 -4 cm is, of course, partially due 
to the fact that the clock is a dynamical experiment and that both components 
of the clock are relatively large compared with, for example, the small rod 
in Beams' experiment; i. e., small components are not necessary. 

Another potential source of error is due to the mode of operation — free 
active axis for the clock as opposed to null centering in the Eotvos experiment. 
As a result, the capacitance controller will tend to produce cross-talk between 
the passive axes and the active axis of the system when the sphere moves off 
the equilibrium position along the active axis. These forces arise when the 
sphere moves out of its equilibrium position and hence its symmetry position 
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with respect to the capacitors, thus giving rise to a restoring force. However, 
this type of force, although much larger than the gravitational force, can be 
compensated for electronically at the same level of measurement accuracy of 

the controller itself and, therefore, will not affect the gravitational restoring 

force and its clock mode of operation. 

The main source of error in the experiment comes not from the deter- 
mination oHhe-period of-oscillation but from the accurate determination of f 

c 

given in equation (21) . If we assume the conservative estimates of others, all 
quantities such as mass, density, and dimension can be obtained to accuracies 
of about 1 part in 10 8 . 


IV. CONCLUSIONS 


The flat-plate spherical mass oscillator was considered as a method 
for obtaining G for the following reasons: 

1. It is mathematically the simplest system. 

2. Internal gravitational forces were relatively large, 

3. The control and measurement systems are relatively simple because 
one member of the oscillator is spherical. 

To this list we could add that the dimensions of the flat plate can als£Lhs_ 
fabricated more accurately than for spheres. Thus, for a given accuracy of 
measuring lengths, one has an inherently more accurate determination of G. 

The effects of density variations can 1)8 me asured and be ultimately compen- 
sated for in calculations of G. 

Since the experiment is relatively simple, the experiment can be dup- 
licated with different substances (change p) , and different masses (change c, 
b, b 1 , t). One then obtains parametric curves co 2 (G, p, c, b, b 1 , r) from 
which G can be extracted with a higher degree of confidence than from any 
individual experiment. 

When one has confidence in a more accurately known G, one can per- 
form Kreuzer-type experiments [ 3 1 3 that can experimentally determine the 
equivalence of active and passive gravitational mass by simply constructing the 
two parts of the oscillator out of different substances. (Compare this with the 
Eotvos experiment that compares inertial versus gravitational mass. ) 
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Finally, even without using experimental techniques comparable to 
those necessary in the satellite superconducting gyroscope experiment or that 
proposed in the Eotvos experiment of Worden and Everitt, accuracies of 1 part 
in 10 6 would easily be obtained without concern for more than fabrication 
accuracies and measurement of the period of oscillation. However, measuring 
density fluctuations and correcting for their effects on the periods of oscilla- 
tions, taking the orbital eccentricity into account, using large components, low 
temperature, and superconductivity, and integrating observations over long 
periods by use of stable atomic clocks would increase the accuracy conserva- 
tively by at least two or three orders of magnitude. 

At the level of accuracy of 1 part in 10 8 or better, it then becomes 
important to consider spatial variations inG or failures in inverse square law. 
To disentangle these effects may require the use of orbits of high eccentricity, 
solar satellites, or both. At any rate, a second independent experimental 
determination of variations of G, presumably by accurately determining the 
magnitude of G as a function of both space and time, is of fundamental 
importance to science. 
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